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Dependence  of  the  Induced  Loss  Factor  on  the  Coupling  Forms 
and  Coupling  Strengths:  Energy  Analysis 


G.  Maidanik 
NSWCCD  (DTMB) 

9500  MacArthur  Blvd. 
WestBethesda,MD  20817-5700,  USA 


Abstract 


In  two  recent  papers  the  induced  loss  factor  is  determined  via  the  modification  to  the 
loss  factor  in  the  linear  impedance  of  a  master  oscillator  caused  by  its  coupling  to  a 
set  of  satellite  oscillators.  A  recapitulation  of  this  linear  impedance  analysis  (LIA)  is 
presented.  A  loss  factor  is  basically  an  energetic  quantity  and,  therefore,  one  may 
inquire  whether  the  induced  loss  factor  may  be  estimated  via  an  energy  analysis 
(EA).  An  answer  to  this  question  is  sought.  It  is  shown  that  the  two  analyses;  (LIA) 
and  (EA),  yield,  in  the  appropriate  frequency  range,  identical  results  for  the  induced 
loss  factor.  This  frequency  range  spans  the  distribution  of  resonance  frequencies  of 
the  satellite  oscillators.  In  this  frequency  range  the  identity  of  the  results  is  not  only 
in  terms  of  gross  features,  but  also  in  details.  Finally,  the  relationship  of  (EA)  to 
statistical  energy  analysis  (SEA)  is  explored.  The  loss  factors  of  the  satellite 
oscillators  are  directly  related  to  modal  overlap  parameters.  It  is  found  that  for  the 
stronger  couplings,  the  validation  of  (SEA)  may  require  these  modal  overlap  para¬ 
meters  to  exceed  a  certain  threshold. 


L  Introduction 


In  two  recent  papers  the  influence  on  the  response  behavior  of  a  master 
oscillator  due  to  its  coupling  to  a  set  of  satellite  oscillators  is  derived  and  examined 
[1,2].  The  complex  is  sketched  in  Fig.  1.  The  examination  is  conducted  in  terms  of 
the  induced  loss  factor  7^(7),  where  (y)  is  the  normalized  frequency  variable; 
y={(ola)Q).  The  normalizing  frequency  is  the  resonance  frequency  (co^)  of  the 
master  oscillator  in  isolation.  The  induced  loss  factor  7^(3;)  accounts  for  the 
modification  in  the  loss  factor  in  the  linear  impedance  of  the  master  oscillator  caused 
by  its  coupling  to  that  set  of  satellite  oscillators.  This  modification  is  the  focus  of  the 
investigation  performed  in  References  1  and  2  and  is  also  the  focus  of  the 
investigation  in  this  paper.  A  recapitulation  of  the  linear  impedance  analysis  (LIA) 
covered  in  these  references  may,  thus,  be  in  order. 

As  Fig.  1  indicates,  the  coupling  between  the  master  oscillator  and  a  satellite 
oscillator  allows  for  stiffiiess,  mass,  gyroscopic  and  a  mixture  of  these  forms  [3, 4]. 
With  the  assistance  of  Fig.  1,  the  linear  equations  of  motion  for  a  master  oscillator  in 
isolation,  for  a  coupled  master  oscillator  and  for  a  typical  satellite  oscillator  are  stated 


in  the  forms: 
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Z»F°(J)  =  75,(3;), 

(la) 

Z„(3')  VoW  =  Pe(y)^ 

(lb) 

r,(y)  =  B,(,y)  V^{y), 

(Ic) 

respectively,  where 

Z"(y)  =  (io>M)li-(yr^  1 

(2a) 

ZJy)  =  (mM)[l-(yr^  {[1  -  5(y)]+  /[% +%(;')]}] . 

(2b) 

R 

1 

*  [^cr  ~  ^cr  )  ]  ~  ^  } 

•  [(1  +  ^cr  )  -  i^iT  f  (1  +  i Vrr  )]"'  }  ^  (2c) 

Br  ( j)  =  -  Wcr  +  i^cr  f  ^ Icr  ) “  ^ iSr  >  j)  ]  [(^  +  ^cr  )  “ (^rr  f  ^ Irr  )  ]~^  (2d) 

(9cr  /  7)^  =  4  m^r  i^cr  f  (1  +  ^  ^cr  )  +  (^r  ^  yf  ’  (2®) 

f  (1 + i  TJ^  )  =  iZrf{\+i7]^)+ f  (1  +  i  7]^r ) .  (2f) 

(z^^  f  =  / yf ;  (z^  =  (x^  /  ;  (z^^  =  (x^^  /  ,  (2g) 

(X^)^  (1  +  /;?^)  =  i.K  i^crf  (l  +  ^>7cr)  =  iKr  >  ^r)  ’  (2h) 

mr  =  {m^/M)‘,  mcr  =  {mcrlm^)-,  gr  =  [Gr /{(Oom^)].  (2i) 
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In  Eq.  (1)  Z^(y)  and  Vf(y),  and  Z^(y)  and  V^iy)  are  the  unmodified  and  the 
modified  linear  impedances  and  responses  of  the  master  oscillator,  respectively, 
Vr(y)  is  the  response  of  the  (r)th  satellite  oscillator  and  Pg(y)  is  the  external  drive; 
this  external  drive  is  applied  to  the  master  oscillator  only  [1,2].  The  other  quantities 
and  parameters  are  defined  in  Fig.  1  and  in  Eq.  (2).  It  is  convenient,  with  only  a 
slight  loss  in  generality,  to  impose  a  similarity  on  the  stiffhesses  (springs)  that  are 
placed  on  either  side  of  the  mass  (m^ )  of  the  (r)th  satellite  oscillator;  namely 

(x,  f  =:  a,  (x^  f ;  {x,,  f  =  (x^ )^  (3) 

where  (x^)  defines  a  designed  distribution  for  the  normalized  resonance  frequencies 
of  the  satellite  oscillators.  In  that  design  the  resonance  fi-equencies  ascend  according 

to  the  value  of  the  index  (r);  i.e.,  (x^^)^  <  (x^qf' ;  q  =  (r  +  1);  1  <  r  <  (i?-l),  and 
the  numbers  of  resonance  frequencies  on  either  side  of  the  resonance  frequency  (<y^) 

of  the  master  oscillator  in  isolation,  are  equal.  As  in  References  1  and  2,  (x^)is 
assigned  the  specific  form 

(:c“)  =  [l+{(l-2F)K«))  r{RMrK2(R)l  (4a) 

where 

F=r(R  +  l)"^  R=R{R  +  \y^\  y=0.6 .  (4b) 

Other  forms  for  (x^)  are  not  only  permissible,  they  are  readily  introduced.  Finally,  it 
is  recognized  that  the  coupling  of  the  (r)th  satellite  oscillator  to  the  master  oscillator 
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is  specified  in  terms  of  the  stiffiiess  parameter  the  gyroscopic  coupling 

parameter  (f^)  and  the  mass  coupling  parameter  {fh^r)-  Th®  values  of  these 
parameters,  either  individually  or  in  combinations,  determine  the  coupling  strengths; 
there  is,  by  definition,  no  coupling  between  a  satellite  oscillator  and  another. 
[cfFig.l.] 

The  induced  loss  factor  ri^{y) ,  investigated  via  Eqs.  (2b)  and  (2c),  reveals  the 

dependence  of  this  quantity  on  parameters  that  define  the  externally  driven  complex; 
see  Fig.  1 .  The  dependence  of  this  induced  loss  factor  on  the  coupling  forms  and 
coupling  strengths  is  of  particular  interest  in  this  investigation.  The  dependence  of 
this  quantity  on  the  distribution  of  the  resonance  firequencies  and  on  the  values  of  the 
modal  overlap  parameters  that  are  assigned  to  the  satellite  oscillators  is  also  of 
interest  in  this  investigation.  Computational  data  of  these  dependencies  are  depicted 
in  a  number  of  figures  in  which  is  presented,  as  a  function  of  (y).  In 

particular,  in  these  figures,  the  phenomenon  of  erosion  is  defined  and  demonstrated. 
In  this  demonstration  the  corresponding  first  order  approximations  to  77^  (y)  play  a 
role  [1,  2].  (The  first  order  approximation  is  evaluated  by  replacing  the  summation 
in  Eq.  (2c)  by  an  integration  and  carrying  out  the  integration  subjected  to  vanishing 
values  of  the  modal  overlap  parameters.)  In  the  evaluations  presented,  certain 
parametric  simplifications  are  introduced.  Largely,  these  simplifications  manifest 
dropping  the  dependence  of  some  parameters  on  the  index  (r);  namely 
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m^=m\  a^^a\  g^  =  g',  W=K  =  b-  (5a) 

Moreover,  to  ensure  that  the  distribution  of  the  resonance  frequencies  span  a 

definitive  and  a  constant  range,  as  a  function  of  (F)  [F=r  (/?  + 1)“’.],  the  parametric 
equality 

(l  +  iw^)  =  (a+a^).  (5b) 

is  conditioned.  Then,  under  the  similarity  conditions  specified  in  Eq.  (3)  and  the 
simplifications  and  condition  introduced  in  Eqs.  (5a)  and  (5b),  respectively,  the  loss 
factors  assigned  to  the  (r)th  satellite  oscillator  are  equal;  namely 

nrr=’Jr  =  ’!cr  =  '!<ry,  ];  (5c) 

where  (x^)  and  y(R)  are  stated  in  Eq.  (4)  [1,2].  Although  ?]{F)  is  not  a  function  of 
the  normalized  frequency  {y) ,  simations  arise  in  which  notationally  it  is  conducive 
to  cast 

n(r)-nr(y)  =  (bimy?r(R)(4fl  4-=<.4iy)-  (Sd) 

The  loss  factor  //(F)  [or  ;7^(>^)],  as  a  function  of  (F),  for  three  values  of  (6); 
6 =(0.1),  (2.0)  and  (10),  is  depicted  in  Fig.  2.  With  the  impositions  stated  in  Eq.  (5) 
the  first  order  approximation  (FOA)  of  the  induced  loss  factor  rj^{y)  assumes  the 
simple  form 

rj,{y)  =  D[C+0{i](,yy\,  7^,(y)^DC,  (6) 

where 

D  =  !  [2r(R))][M ,  /(MR)  [l+m,T' ;  M,  =  (Rm) 


(7a) 
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C  =  m,+a,fHglyf],  (7b) 

(9  =  [(l  +  m^-ac)«c]»  ('7c) 

and  the  term  0{ri{yyf'\s  suppressed  and  lumped  with  the  higher  order  approx¬ 
imations  [2].  The  validity  of  Eq.  (6)  is  predicated  on  the  condition  that 

’l,^y)  =  n<7o)  =  WR  +  \T'][(yfr(R)]r=r,(.yy,  z”=l; 

[i+(r/2)r<''^>  <j<[i-(r/2)p''7).  (8) 

[cf.  Eq.  5d.)  It  transpires  that  the  first  order  approximation  of  i.e.,  rij{y), 

coincides  with  the  mean-value  averaging  of  the  undulations  that  beset  the  exact 
levels  of  the  induced  loss  factor  rj^iy)  when  the  modal  overlap  parameter  (b)  is  small 
compared  with  unity  [1,  2,  5].  [cf  Fig.  3d  of  Reference  2.]  The  first  order 
approximation  rii{y)  of  the  induced  loss  factor  is  thus  independent  of  (b).  Then,  the 
first  order  approximation  rij{y)  of  %(y)  is  free  of  erosion.  (It  is  speculated  that  the 
erosion  is  manifested  in  the  higher  order  approximations  to  rj^iy) .  Since  this  higher 
order  approximation  procedure  is  geared  toward  accounting  for  the  higher  values  of 
(b),  the  undulations  at  the  lower  values  of  (b)  are  not  expected  to  be  accounted  for  in 
this  process  of  approximations.)  In  Reference  2  computations  of  T]g{y)dXQ  con¬ 
ducted  and  presented  in  Figs.  3e  and  3f  and  4-7.  Each  figure  in  this  series  pertains 
first  to  a  coupling  form  and  second  to  a  coupling  strength.  Also,  each  figure  is  a 
superimposition  of  four  curves;  three  of  the  four  pertain  to  an  exact  evaluation  of 
ri^{y)  for  three  modal  overlap  parameters;  i.e.,  6  =  (0.1),  (2.0)  and  (10) .  The 
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fourth  curve  pertains  to  the  first  order  approximation  stated  in  Eqs.  (6)  -  (8).  In  this 
manner  the  content  of  information  in  these  figures  is  high  and  can  be  deciphered  at  a 
glance.  It  is  assumed  that  these  figures  are  available  to  the  reader. 

A  loss  factor  is  basically  an  energetic  quantity  and,  therefore,  one  may  inquire 
whether  the  induced  loss  factor  %(y) ,  as  just  determined  via  the  linear  impedance 
analysis  (LIA),  may  be  estimated  via  an  energy  analysis  (EA).  In  this  latter 

determination  the  induced  loss  factor  is  designated  rig{y) .  A  question  arises:  Are 

TJsiy)  and  rf^siy)  5  or  to  that  matter  the  respective  first  order  approximations  rjj  (y) 

and  ^/(y) ,  respectively,  identical  quantities?  A  sense  of  such  an  identity,  or  even  a 
disparity,  if  any,  may  help  interpret  the  salient  features  that  they  commonly  and 
individually  possess.  To  this  end  the  paper  is  largely  dedicated. 
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II.  Derivation  of  the  Energy  Equation  of  Motion 

The  linear  equation  of  motion  of  the  complex  composed  of  a  master  oscillator 
coupled  to  a  set  of  satellite  oscillators  may  be  largely  derived  via  the  Lagrange 
equations.  The  Lagrangian  describes  the  difference  between  the  kinetic  and  potential 
energies  stored  in  the  oscillators  and  in  the  couplings  [3],  (Notwithstanding  that  the 
linear  equations  are  correctly  stated  in  Reference  3,  a  persistent  typographical  error 
in  the  preceding  Langrange's  equations  needs  to  be  corrected.)  Although  the  kinetic 
and  potential  energies  are  here  determined  separately,  it  is  the  energy  and  not  the 
Langrangian,  as  such,  that  is  of  immediate  interest.  The  kinetie  energy  (y )  plus 
the  potential  energy  E^piy)  stored  in  the  master  oscillator  is  given  by 

EQ{y)=Egf^(^y)  +  E^piy)',  E^f^^y)  =  {\I2)M  \VQ{y)\  , 

Eop (y)  =  (yr^  (t) ;  y=ioj/co^),  (9) 

where  E^^iy)  is  the  stored  energy  in  the  master  oscillator  [3].  Similarly,  the  kinetic 
energy  E^^iy)  ^h®  potential  energy  that  is  stored  in  the  (r)th  satellite 


oscillator  is  given  by 

E,  (y)  =  E,k  (y)  +  E,p  (y) ;  E,k  (y) = (1  /  2)  m J  V,  (y )  p ; 
Efpiy)~i^r^  ^rK^y^’ 


(10) 
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where  Ej.{y)  is  the  stored  energy  in  the  (r)th  satellite  oscillator  [3].  In  addition,  the 
kinetic  energy  ^crK  (y)  and  the  potential  energy  stored  in  the  coupling, 

between  the  master  oscillator  and  the  (r)th  satellite  oscillator,  is  given  by 
EAy)  -  ^crK  (y)  +  ^crP  (y),  (lla) 

EcrKM  =  [m,  |  V,(y}  +  V,(y)  f 

+  Im{  (g,/y)ir^(y)  +  K{y)][V„(y)-V,(,y)i  j]  ,  (lib) 

E„p=(U2)m,  (z,,  f  I  V„(y) -V,(y)f,  (11  c) 

where  E^^  (y)  is  the  stored  energy  in  the  coupling  between  the  master  oscillator  and 
the  (r)th  satellite  oscillator  [3].  Using  Eq.  (2)  one  obtains,  from  Eqs.  (9)  -  (1 1),  the 
expression  for  the  stored  energies 

£„(>>)  s  [1+W"l;  E„(y}^[E„(y)/E,f;(y)l  (12a) 

E,(y)={E,(y)l  E„n{y)]  =  (12b) 

E„(y)^[E„(y)IE„^(yy\ 

=  m,[m„  1 1+B,  +(z„)^  I \-B,  I  y)(B, -B')],  (12c) 

respectively,  where  a  single  bar  over  a  stored  energy  quantity,  e.g.,  Ej.  (y) ,  indicates 
a  normalization  by  the  energy  quantity  EoK  (y) .  In  particular,  the  so  normalized 

energy  ^'^^(y),  stored  in  the  (r)th  satellite  oscillator  and  in  the  coupling  of  this 
satellite  oscillator  to  the  master  oscillator,  is  given  by 
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_  _  _  _  R  _ 

E„(y)  =  E,(y)  +  E„(y) ;  E„(y)  =  Y.^„,(y),  (US) 

1 

where  Eg^(y)  is  the  normalized  energy  stored  in  the  satellite  oscillators  and  in  the 
couplings  of  these  oscillators  to  the  master  oscillator. 
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III.  Conservation  of  Energy  -  the  Balance  of  Power  - 
and  the  Derivation  of  the  Induced  Loss  Factor 

Remembering  that  the  external  drive  acts,  by  definition,  only  on  the  master 
oscillator,  the  normalized  external  input  power  II ^iy)  that  is  available  to  maintain 
the  stored  energies  in  the  complex  may  be  derived  from  the  linear  equation  of 
motion.  This  equation  is 

Z<,(®)F„(ffl)  =  />,(«),  (13a) 

,  (13b) 

where  Z^(ru),  Vq{o))  and  P^{o))  are,  respectively,  the  linear  impedance,  the 
response  and  the  external  drive  that  induces  the  response  in  the  master  oscillator 
when  coupled,  [cf.  Eq.  (lb)  and  (2b)  and,  also  Eqs.  (5)  and  (7a)  of  Reference  2.] 
This  external  input  power  is  derived  in  the  form 

n,(y)  =  [Tl,(y)l{<oE„^(,y))]  -  n,(,y)  =  Re{/>,(y) 

(y^/2)njy)  =  [^„  +  T!,(y)].l,  (14) 

where,  hereafter,  a  single  bar  over  a  power  quantity;  e.g.,  II ^iy),  indicates  a 
normalization  by  the  power  quantity  {(x>EQj^{y)}.  [cf  Eq.  (12).]  The  third  of  Eq. 
(14)  states  that  a  portion  no(>’)  of  H^iy)  is  dissipated  in  the  master  oscillator  such 
that 


(y^  I2)n„(y)  =  nl{y)-y,  <lHy)  =  no-.  n<,(;v)=[n„(;-)/{®£„^(y)}],  (15a) 
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and  the  portion  noj(>’)  of  H^iy)  is  dissipated  in  the  satellite  oscillators  and  in  the 
couplings 

(y^l2)n^(y)=4W-U  n„,(y)  =  [n„,(y)l{coE„^{,y))]  (15b) 

where,  as  already  intimated,  the  superscript  (e)  indicates  that  the  loss  factors  are 
determined  via  an  energy  analysis  (EA)  and  not  via  a  linear  impedance  analysis 
(LIA).  The  conservation  of  energy  (the  balance  of  power)  demands  that 

no(y)+njy)  =  n,(y).  (16) 

Employing  Eq.  (12d),  Eq.  (15b)  may  be  decomposed  in  the  form 

_  _  _  _  _ 

n^(y)=Yn„iyy,  n„(y)  =  n,(y)  +  ii,Ay).  ( 1 7) 

1 

where  by  definition 

(y )  =  T],E^  (y) ;  (y)  =  17^,  (j)  .  (1 8) 

The  loss  factors  (77^)  and  (77^^)  are  the  stiffness  control  loss  factors  associated  with 
the  springs  on  the  fore  and  the  back  sides  of  the  mass  (777^)  of  the  (r)th  satellite 
oscillator.  The  spring  on  the  back  side  constitutes  the  stiffiiess  coupling  form.  The 
spring  on  the  fore  side  renders  the  satellite  oscillator  an  oscillator  rather  than  merely 
a  mass.  It  is  assumed  that  the  damping  is  provided  in  these  springs  only  [1,  2]. 
Provisions  for  other  types  of  damping  can  be  made,  but  the  increase  in  algebraic 
complexity  can  hardly  be  justified  at  this  stage. 
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From  Eqs.  (15b),  (17)  and  (18)  one  obtains  the  energetic  version  of  the 

induced  loss  factor  ijI  (7)  in  the  form 
R  _ 

rjl(y)  =  (/  E,(y)  +  7,,  E„(y)  )  ,  (19) 

I 

where  Ej.{y)  and  E^.^{y)  are  more  explicitly  expressed  in  Eqs.  (12b)  and  (12c).  The 
induced  loss  factor  (y) ,  derived  via  (LIA),  is  explicitly  expressed  in  Eqs.  (2b)  and 

(2c)  [1,2].  To  what  degree  then  is  Tj^(y) ,  stated  in  Eq.  (19),  equal  to  t]^  (y)  ?  To 
establish  analytically  the  answer  to  this  question  may  require  undue  algebraic 
manipulations  which,  again,  can  hardly  be  justified  at  this  stage.  Instead,  the 
equivalence  is  cursorily  tested  computationally.  Imposing  the  simplifications  stated 

in  Eq.  (5),  the  explicit  expression  for  tj^ (y)  may  be  reduced  to  the  form 
R  _  _ 

r,:(y)  =  (y^/2)J^^iF)E„(y)^(y^/2)n‘(y)E„(y)  ,  (20) 

1 

where  E^^iy)  is  stated  in  Eq.  (12d),  Tj^(y)  is  the  loss  factor,  as  a  function  of  (y),  of 
a  typical  satellite  oscillator  with  its  resonance  frequency  in  the  vicinity  of  (y) ,  tj(F) 
is  stated  in  Eq.  (5c)  and  depicted  in  Fig.  2,  and  E^^iy)  is  the  normalized  energy 
stored  in  the  (r)th  satellite  oscillator  and  in  the  coupling  of  this  oscillator  to  the 


master  oscillator;  namely 
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E,,iy)  =  in,  (l  +  m,r2{  f  +  W)f  }"* 

•  {{\+a{z^f  ][{m^-\-a^{z°ff  +  {a^T]{r){z°f  -  {g,  !  y)f] 

+  [«c  (^r  ]  [{(l+2m,)-a  (z^f  f  +  {L(y)f  ] 

-(gr^y)[U-(i  +  ^c+^c)(^rf}(^^(^)(^rf  +  (Sr  I  y)) 

-  {{\+2m^)-a{z^,f)  {{Um^+ac)  ri{r)-{g, !  y)}]} ,  (21a) 

where  the  explicit  expressions  for  ?+(>')  and  t_{y)  are 
t+iy)  =  (1  +  Wc)^(^)(^r)^+ ; 

/_  ( j)  =  (1  +  )^  -  [  <^c  (^r  )^  -  (^r  /  J)]  •  (21b) 

Specific  parameters  involved  in  Eqs.  (20)  and  (21)  are  defined  in  Eqs.  (2)  —  (5). 
Using  Eqs.  (20)  and  (21)and  the  simplifications  and  condition  stated  in  Eq.  (5), 

rj^iy)  is  computed  and  the  results  are  presented  in  Figs.  3-6.  These  figures  are  cast 
in  the  format  of  Figs.  3e  and  3f  and  4-7  of  Reference  2.  Thus,  each  figure  is  the 

superimposition  of  four  curves:  three  exact  evaluations  of  rj^iy)  for  the  three  values 
of  the  modal  overlap  parameter  (b);  b  =  (0.1),  (2.0)  and  (10)  and  one  for  the  first 

order  approximation.  (The  first  order  approximation  for  rj^iy) ,  that  is  employed,  is 
derived  subsequently.)  Again,  in  the  manner  of  this  display  the  content  of 
information  in  these  figures  is  high  and  can  be  deciphered  at  a  glance. 
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notwithstanding  that  the  comparison  with  Figs.  3e  and  3f  and  4-7  of  Reference  2  is, 
thereby,  facilitated. 

In  References  1  and  2  the  summation  for  (y)  is  replaced  by  an  integration. 
The  first  order  approximation  of  this  integral  is  evaluated  and  is  stated  in  Eqs. 
(6)  -  (8).  A  similar  procedure,  replacing  the  summation  by  an  integration,  is  now 
applied  to  Eq.  (20).  The  first  order  approximation  only  is  evaluated.  The  induced 
loss  factor  so  evaluated  yields 

r,‘,(y)  =  DC,  (22) 

where  (D),  (C)  and  the  range  of  validity  in  (y)  are  as  stated  in  Eqs.  (7)  and  (8),  but 

0^=(l/2)[(l-i-mc)(^c+«c)+2w^ac]-  (23) 

In  the  frequency  range  of  validity,  the  loss  factors  ri{y)  and  r}^{y)  in  Eqs.  (6)  and 
(22),  respectively,  are  found  to  be  identical.  From  Eq.  (8)  one  finds 

'iiy)  =  ’i"W={b(Ry^r'][(yfriR)V, 

[  1  +  (y  /  2)  r*' '  S  >-  S  ( 1  -  (r  /  2)  .  (24) 

The  evaluations  of  T]^{y)  [and  T]{y)'\  for  b  =  (0.1),  (2.0)  and  (10)  are  presented 
graphically  in  Fig.  7.  In  Fig.  7a,  R=27  and  in  Fig.  7b,  R=7.  The  essential  similarity 

of  7]^  iy)  with  the  corresponding  t](F)  is  clear  when  Figs.  7  and  2  are  appropriately 
compared.  (One  recalls  that  the  loss  factor,  presented  graphically  in  Fig.  2,  is  the 
loss  factor  that  is  designed  and  assigned  to  the  (r)the  satellite  oscillator  [2].)  Since 

ri{y)  and  T]^{y)  are  identical,  the  first  order  approximation  (FOA)  of  ri^{y)  and  of 
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Tis(y),  stated  in  Eqs.  (6)  and  (22),  respectively,  differ  in  that  the  coefficients  (O) 

and  (O^)  are  not  identical.  Usually,  the  terms  (DO)  {;7(y)}^and  (DO^)  {rj^(y)}^ 
are  negligible  and,  therefore,  the  difference  between  them  is  rarely  significant. 
Notwithstanding  that  these  terms  are  considered  to  belong  with  the  higher  order 
approximations  and,  as  such,  they  are  removed  from  the  first  order  approximations 

[2].  Thus,  it  is  concluded  that  the  first  order  approximations  of and  //JCy)  are 
identical.  Moreover,  from  Eqs.  (6)  and  (22)  these  first  order  approximations  are 
independent  of  the  modal  overlap  parameter  (b)  [2]. 

Remarkably,  the  equivalence  between  Figs.  3-6  and  the  corresponding 
Figs.  3e  and  3f  and  4-7  of  Reference  2  is  not  only  in  terms  of  gross  features,  but  also 
in  details;  e.g.,  the  erosions  in  Figs.  3e  and  3f  and  4-7  of  referenee  are  duplicated,  to 
a  tee,  in  Figs.  3-6,  respectively,  at  least,  in  the  range  of  the  normalized  frequency  of 
concern,  [cf.  Eqs.  (8)  and  (24).] 
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rV.  Renormalization  of  Stored  Energies  and  Powers  - 
Illegitimate  and  Legitimate  Loss  Factors 

It  has  been  adopted  that  quantities;  e.g.,  the  stored  energy  and  the 

dissipated  power  11^^  (y),  when  normalized  by  ^oK  (>')and  {a>  EoK(y)). 

respectively,  are  designated  by  a  single  bar;  namely,  and  110^(7) .  Situations 

arise  in  which  a  normalization  by  E^iy)  and  {(oE^iy)},  respectively,  may  be 
preferred.  Such  normalizations  are  to  be  designated  by  a  double  bar,  namely, 

^03(7)  and  n^^(y),  respectively.  In  particular  then 

nM  =  [n„W/!®£„W}]  =  n„w  =  2[i+(yy‘ )"'  >jo,  (zsa) 

nos(y)  =  [Ti„{y)i{coE,{y)}\  =  =  2[iHyf  r‘  %‘W,  (25b) 

nlsiy)  =  'E,’7(E)Eor(y)  =  'l"MS„(y).  {25c) 

1 

It  emerges  then  that  whereas  the  loss  factors  (7^,)  and  (77J)  relate  to  the  normal¬ 
ization  employing  the  kinetic  energy  E^^^iy)  and  the  power  {oE^fr  (y)},  the  loss 

factors  //ooCy)  and  r]lsiy)  relate  to  the  corresponding  quantities  normalized  by  the 
stored  energy  E^iy)  and  the  power  {coEQ{y)).  The  stored  energy  is  the 

combined  energy,  consisting  of  the  kinetic  energy  EoK  (y)  and  the  potential  energy 
Eopiy)  stored  in  the  master  oscillator;  E^^y)  =  EqK  (y)  +  £^p(y),  as  stated  in 
Eq.  (9).  Equation  (25c)  is  then  a  renormalized  version  of  Eq.  (20).  The  ratio  of  the 
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energy  (y)  stored  in  the  satellite  oscillators  and  in  the  couplings  to  that  of  the 
energy  E^iy)  stored  in  the  master  oscillator  may  be  dubbed  the  global  coupling 


strength  and,  as  such,  is  designated  .  From  Eqs.  (25a)  and  (25b)  one  obtains 


'7«W=2  [\Hyf  ]  '  vtW  =  Wiy)  sJW] ; 


K(y)  =  Eos(y)  =  [E„(y)l  E,(.y)]. 

Were  one  to  construct,  from  Eqs.  (22)  and  (26a)  the  hybrid  expression 


[;7*(y)  =*(>.)]*  =  2[i+yr' 


b<\, 

ysiy)’  b>\. 


(26a) 


(26b) 

(26c) 


then,  except  for  erosion  at  the  higher  values  of  the  modal  overlap  parameter;  b>\, 
one  finds  the  hybrid  quantity  to  be  largely  undulations  free  and  independent  of 

(b).  Referring  to  Eqs.  (26b)  and  (26c),  it  follows  that  in  the  hybrid  milieu,  7^(y)  is 

inversely  proportional  to  Ho(y)  and  vice  versa.  Taking  note  of  Eq.  (24),  there  is  no 

way  for  a  physically  acceptable  complex  to  entertain  a  vanishing  rj^iy) .  Of  course, 
Eq.  (26a)  does  not  fall  under  this  spell,  but  neither  can  one  cavalierly  ignore  the 

undulations  that  beset  (y)  when  (b)  is  reduced  with  the  intention  of  rendering 


7^(y)  negligible  [2]. 
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Being  cognizant  of  the  second  of  Eq.  (26a),  the  global  coupling  strength 
may  feature  in  another  renormalization  in  which  the  entire  energy  Eg{y) 


stored  in  the  complex  is  employed.  This  stored  energy  is  given  by 

Eeiy)  =  E,{y)^E,,{y)  -  E,{y)  [l+Kiy)],  (27) 

and  the  corresponding  renormalization  for  the  powers  is  given  by  {®£'g(>^)} .  Such 
renormalizations  are  to  be  designated  by  a  tilde  over  the  quantity;  e.g. 

Ksiy)  <E,,{y)IE,{y)\,  U,,{y)={Yl,,{y)l{(DE,{y)}\  (28a) 

From  Eqs.  (25),  (27),  and  (28a)  one  may  derive 

Eos(y)  =[£„(:i')/£.W]=f„W[l+H*(y)r',  (28b) 

n„(y)  =  [n„,(y)l{0lE,(yy,]  =  n„,(y)[l  +  H'(j.)r' ,  (28c) 

!,(;.)  =  [1+Hj(y)].  (28d) 


The  two  renormalizations  just  introduced,  in  Eqs.  (25)  and  (28),  lead  to  the 
definition  of  two  distinct  loss  factors.  The  first  is  defined  in  terms  of  the  energy 
Eoiy)  stored  in  the  master  oscillator  only;  namely 

(t)  =  {^e{y)l{o)E^{y)}]  =  [7oo(t)+^^(t)  So(t)].  (29a) 

where  use  is  made  of  Eqs.  (14)  and  (16)  .  The  second  is  defined  in  terms  of  the 
energy  E^{y)  stored  in  the  entire  complex;  namely 

n‘Ay)  =  U,(y)l{a,E,(y)}  =  ^fWll+s;)^)]-'.  {30a) 
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where  use  is  made  of  Eqs.  (28)  and  (29a).  These  two  loss  factors  were  previously 

defined  and  discussed  [6,  7].  In  these  discussions  it  was  claimed  that  is  not  a 
legitimate  loss  factor.  This  designation  was  primarily  predicated  on  assigning  to  the 
definition  of  this  loss  factor  the  entire  external  input  power  Hgiy) ,  but  accounting 
for  the  stored  energy  ^^^(y)  that  is  maintained  in  the  master  oscillator  only,  thereby, 
ignoring  the  portion  of  the  stored  energy  that  is  maintained  in  the  satellite  oscillators 
and  in  the  couplings.  (No  wonder  the  question  of  "where  did  the  energy  go?"  arose 

from  accepting  77f(>')as  a  legitimate  loss  factor  [8-13].)  In  this  connection  one 
recalls  that 

vfiy)  =  [7oo(T)+7o.(T)];  vlsiy)  =  [ri\y)Kiyn  (29b) 

Again,  since  7]%{y)  is  found  to  be  largely  independent  of  Tj^{y),  one  concludes, 
from  Eq.  (29b),  that  Tjf(y)  is  also  independent  of  Tj^(y) .  Equation  (29)  makes  clear 

that  to  change  ?jf(y)  either  the  coupling  parameters,  the  mass  ratio  or 

both,  need  changing.  In  this  assessment,  is  assumed  to  be  fixed,  and 

7oo(T)  =  27o(l+y^)~*-  [cf.  Eq.  (25a).]  On  the  other  hand,  the  second  loss  factor, 

dubbed  the  effective  loss  factor  and  designated  rj^iy),  does  take  into  account  the 
whole  energy  stored  in  the  complex  [6,  7].  It  may  thus  be  designated  a  legitimate 
loss  factor.  From  Eqs.  (27)  and  (30a)  one  finds 
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>?Uy)^'ifM{n‘(y)l’?‘iy)  +  noAy)r'}  s  ifiy)-  (30b) 

In  particular,  if  the  induced  loss  factor  rj^siy)  exceeds  the  indigenous  loss  factor 
Voo(y)  of  tlie  master  oscillator,  Eq.  (30b)  may  be  further  reduced  to 

\nos 

vl ( j)  =  {n\y) (y) ] (3^) +vlsiy)]  <\  ;  (y)  <  ^os (y) •  POc) 

ri\y) 

Thus,  the  effective  loss  factor  ril{y) ,  under  this  condition,  is  a  parallel  combination 

of  the  induced  loss  factor  7o^(y)and  the  loss  factor  rj^iy)  assigned  to  a  typical 

satellite  oscillator.  It  follows  that  ril{y)  is  less  than  either  one  of  these  loss  factors. 

Consider  a  reasonable  complex  for  which  >  ifiy)-  Under  this  additional 

condition,  Eq.  (30c)  may  be  further  reduced  to 

vtiy)=v\y)’  ilsiy)'>n\y')-  (30d) 

Were  a  noise  control  effort  intended  to  achieve  an  effective  loss  factor  7|(y)that 
exceeds  the  indigenous  loss  factor  ri^{y)of  the  basic  master  oscillator,  a  number  of 
design  criteria  need  be  implemented.  In  this  exercise  the  noise  control  is  to  be 
derived  from  the  attachment  of  a  number  of  satellite  oscillators  to  this  master 
oscillator  [12,  15].  [cf.  Fig.  1.]  If  the  two  conditions  stated  in  Eqs.  (30c)  and  (30d) 
are  to  be  satisfied,  it  follows,  from  Eqs.  (22)  and  (27),  that  the  three  design  criteria 

^00  (t)  <  (t)  So (t) = (t)  ,  (31  a) 


(31b) 
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Va,{y)<'n^{y).  (3ic) 

must  be  simultaneously  satisfied  [7].  The  first  of  these  criteria,  then,  demands  that 

the  induced  loss  factor  rf^^iy)  exceeds  the  loss  factor  rf^iy)  of  the  isolated  master 

oscillator.  To  satisfy  this  criterion,  it  is  observed  that  both  rfiy)  and  S^(7)be 
designed  to  possess  high  values,  as  dictated  in  Eq.  (31a).  How^  high?  As  stated  in 
Eq.  (31b),  the  energy  stored  in  the  satellite  oscillators  and  in  the  couplings 

need  be  in  excess  of  the  energy  £'o(y)  stored  in  the  master  oscillator.  This  is 
achieved  by  increasing  the  coupling  parameters  and  the  modal  density  (number  of 
satellite  oscillators  per  unit  frequency)  of  the  satellite  oscillators  [7].  Simul¬ 
taneously,  a  typical  loss  factor  r)^{y)  of  the  satellite  oscillator  need  be  maintained  in 
excess  of  rj^iy) .  This  condition  ensures  that  the  dampings  in  the  satellite  oscillators 
can  handle  the  excess  in  the  stored  energy  that  they  are  designed  to  harbor  [7].  This 
handling  needs  to  exceed  that  ’which  is  indigenous  to  the  master  oscillator,  and  hence 
the  criterion  stated  in  Eq.  (31c).  Details  of  these  criteria  may  be  cast  in  terms  of  the 
parameters  that  describe  the  master  oscillator,  the  satellite  oscillators  and  the 
couplings  between  them.  These  details  may  be  drawn  from  the  preceding  equations. 
To  help  with  the  interpretation  of  these  details,  a  number  of  computations  are 

performed.  In  these  computations  the  loss  factors  rif{y)  and  7]l{y),  as  fimctions  of 

(y),  are  evaluated  and  contrasted.  Again,  only  a  few  representative  cases  are 
displayed.  The  displays  are  given  in  Figs.  8-11.  The  presentation  covers  cases  that 
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parametrically  conform  to  those  governing  Figs.  3-6,  respectively.  The  two  loss 

factors;  rif{y)  and  rj^iy),  are  contrasted  on  the  same  figure  and  each  is  com¬ 
pounded  by  a  superimposition  of  curves  that  pertain,  in  turn,  to  the  modal  overlap 
parameters  b  =  (0.1),  (2.0)  and  (10)).  [cf  Fig.  7.]  In  addition  to  the  evaluations  of 

//f  (y)  and  //Ky),  the  corresponding  levels  of  [(7?)“^  Ho(y)]  are  also  evaluated,  as 
a  function  of  (y) ,  and  are  displayed  in  Figs.  12-15.  Again,  three  curves  that  pertain 
to  the  modal  overlap  parameters  b  =  (0.1),  (2.0)  and  (10)  are  superimposed  in  these 
figures.  The  undulations  in  the  levels  when  (b)  is  less  than  unity  and  the  suppression 
of  the  undulations  when  (b)  is  in  excess  of  unity  is  clearly  apparent  in  these  figures 

too.  The  first  order  approximation  of  [(/?)"’  Ho(y)]  are  also  superimposed  on  these 
figures.  This  superimposition  exposes,  once  again,  the  phenomena  that  the  mean- 
value  averaging  of  the  undulations  coincide  with  the  first  order  approximation  (FOA) 
and  that  erosion  commences  and  increases  as  (b)  approaches  and  increases  beyond 


unity  [1,2,  5,  12]. 


V.  Relationship  to  the  Statistical  Energy  Analysis  (SEA) 


It  may  be  conducive  to  include  in  this  thesis  a  possible  relationship  between 
the  energy  analysis  (EA)  dealt  with  in  the  preceding  three  sections  and  the  statistical 
energy  analysis  (SEA).  The  latter  analysis  was  initiated  in  the  early  1960s  at  BBN 
and  subsequently  has  become  a  major  tool  in  noise  control  engineering  [3,  16].  A 
rudimentary  (SEA)  is  applied  to  a  complex  comprising  of  a  master  oscillator  coupled 
to  a  set  of  individual  satellite  oscillators;  the  satellite  oscillators  are  neither  coupled 
to  each  other  nor  externally  driven,  [cf  Fig.  1 .]  In  terms  of  (SEA)  the  equations  that 
govern  the  energy  distribution  in  the  complex  are 
R  R 

[7oo(t)+Z^/-o(t)]  Eo{y)-Y,VoAy)Eoriy)  =  (32a) 

1  1 

[Vr  (T)  +  lor  (t)  ] ^or  (T) “  Iro  (t)  (t)  =  0  ,  (32b) 

where  E „{}>),  E^^iy),  laiiy)^  Ir^y)  and  ng(y)are  previously  defined  in  Eq.  (9), 
(12d),  (25a),  (5d)  and  (14),  respectively  [3,  16].  In  Eq.  (32)  7o/-(t)  and  ^/^^(y)  are 
the  coupling  loss  factors  from  the  (r)th  satellite  oscillator  to  the  master  oscillator  and 
vice  versa,  respectively.  In  a  conservative  coupling  =  77^^  (y)  [3].  After  a 

straightforward  algebraic  manipulation  of  Eq.  (32)  one  obtains 

R 

it  iy)  =  [?7oo  (y)  +  its  (t)  ];  its  (t)  =  E  (t)  Co  (t)  ; 

1 


Co  (t)  =  [Eor  (t) / Eo (y)]  =  lroiy)ilr  (t)  +  lor  (t)  ]  ’ . 


(33) 
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where  Eq.  (29)  is  consulted  and  ^g(y)is  the  modal  coupling  strength  of  the  (r)th 
satellite  oscillator  to  the  master  oscillator  [3,  16].  From  Eqs.  (29)  and  (33)  one 
further  obtains  that 

R 

nlsW  =  ’i‘(y')K(.y)<  (34a) 

1 

In  this  role  the  stored  energy  ratio  E^iy)  may  be  designated,  as  already  mentioned, 
the  global  coupling  strength.  Again,  with  a  slight  of  statistical  hand,  Eq.  (34a)  is 
approximated 

Kiy)  =  C^y)  {N,{y)U{R-N,iy)}  +  i?t/{Ar,(3;)-;?}],  (34b) 

where  Ng{y)  is  the  number  of  satellite  oscillators  that  contribute  viable  modal  coupling 

strengths.  A  typical  modal  coupling  strength  is  designated  (y) ;  the  viable  modal 
coupling  strengths  are  assumed  to  be  statistically  indistinguishable;  The  use  of 
viability  is  in  reference  to  investigations  revealing  that  the  coupling  loss  factor, 
typically  rio^(y)  or  JJro(y)^  decreases  as  the  frequency  disparity  increases  between 
the  resonance  frequency  of  a  satellite  oscillator  and  the  resonance  frequency  of  the 
master  oscillator  to  which  it  is  coupled  [3].  As  the  disparity  between  the  resonance 
frequencies  increases  beyond  a  specific  threshold,  the  modal  coupling  strength 
decreases  sharply  [3,  4].  The  threshold  is  set  by  either  a  critical  value  of  the 
frequency  disparity  when  typically  r}^(y)  >  rj^^iy)  or  a  critical  value  of  the 
frequency  of  transition  (y^)  when  typically  ?j^{y)  <  Tj^^{y), both  inequalities  are 
validated  only  in  a  frequency  range  for  which  (y)  is  at  and  in  the  vicinity  of  unity.  In 
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the  latter  case  the  frequency  of  transition  is  defined  by  (y^)  = (y,) . 
Consulting  Reference  1  and  imposing  the  simplification  rendered  in  Eq.  (5a),  one 
may  derive 

Ns{y)=^{^0)b)[(D,ri\y)Y^  ==(i?  +  l)(Ay/y)[(l  +  w^)*^2(y2)y(^)]-*.  (35a) 

Using  Eqs.  (4)  and  (24)  one  finds  that 

Y{R)  =  {yl{2R)]\  (Ay/y)  =  (y/2);  y  <  1,  (36a) 

and,  therefore,  from  Eqs  (35a)  one  obtains 

N,(y)  =  R  [(l  +  m,)‘'2  (y)^r'  =  R,  (35b) 

indicating  that  all  the  satellite  oscillators  which  resonance  frequency  distribution  is 
defined  in  Eqs.  (2)  and  (4),  may  be  designated  viable.  Equations  (34b)  and  (35b) 
yield 

C(y)  w' K(yy.  S^Cj.)  =  (^)  (37) 

Thus,  the  typical  modal  coupling  strength  ^o(y)  is  one  (l/i?)ththe  global  coupling 
strength,  both  quantities  are  appropriately  averaged  [3].  From  Eqs.  (20),  (25),  (26) 
and  (37)  one  obtains 

C(>')  =  (y)^[(l  +  7'^)«r'i;£„.W,  (38a) 

1 

where  Eo^(y)is  explicitly  stated  in  Eq.  (21).  Under  the  condition  that  allows  the 
summation  in  Eq.  (19)  to  be  replaced  by  integration,  one  derives  for  the  first  order 
approximation  of  a  typical  modal  coupling  strength 
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Ciy)  =>(l(y)  =  (4'!-/*)  (M,IM)[r^/(l  +  y^)  (l  +  m,)]-'(C); 

C  =  [(m^+aj)^  +(?/y)^); 

(l  +  (y/2)r^‘'^>  <  y  S  [l-(r/2)r<‘'^> ,  (38b) 

where  use  is  made  of  Eqs.  (6),  (7),  (22)  and  (23).  A  number  of  sets  of  parametric 
values  that  define  previous  figures  are  identified  and  applied  to  Eq.  (38a).  The 
results  of  this  application  are  presented,  for  three  values  of  (b) ;  namely,  b  =  (0. 1), 
(2.0)  and  (10),  in  Figs.  12-15.  Also  superimposed  in  these  figures  are  the  levels  of 

the  first  order  approximation  of  Co(y)  ’  stated  in  Eq.  (38b).  Moreover,  if  there  are 
undulations  in  levels  in  curves  based  on  Eq.  (38a),  these  undulations  are  suppressed 
in  the  corresponding  curves  based  on  Eq.  (38b).  Again,  neglecting  the  phenomenon 
of  erosion,  this  suppression  is  commensurate  with  the  appropriate  mean-value 
averaging,  notwithstanding  that  the  phenomenon  of  erosion  is  found  in  these  figures 
too.  An  immediate  violation,  however,  is  revealed  in  Figs.  12-15;  in  particular,  in 
figures  pertaining  to  strong  couplings  and  a  modal  overlap  parameters  (^)that  are 

small.  Since  the  maximum  value  of  a  modal  coupling  strength  ^^(y)in  a  rudi¬ 
mentary  (SEA);  e.g.,  in  Eq.  (33),  cannot,  by  definition,  exceed  unity,  (SEA)  is 

invalidated  for  a  complex  for  which  a  typical  modal  coupling  strength  (y) ,  as 
stated  in  Eq.  (38),  exceeds  unity.  For  example,  the  levels  pertaining  io  b  =  (0.1)  in 

Figs.  12a,  13a  and  15a,  which  display  (^o(y) ,  as  a  function  of  (y),  stand  in  obvious 
violation  of  a  core  tenet  of  (SEA).  The  contradictions  that  arise  in  Fig.  12a  and 
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others,  in  which  levels  of  <^0(7)  that  are  in  excess  of  unity  are  found,  do  not, 
however,  negate  the  analysis  under  which  they  are  derived,  they  merely  lie  outside 
the  range  of  validity  for  the  use  of  (SEA).  It,  thus,  emerges  that  (SEA)  has  a  modal 
overlap  parameter  threshold;  situations  arise  in  which  (SEA)  is  not  valid  for  certain 
degrees  of  coupling  strengths  unless  the  modal  overlap  parameters  assigned  to  the 
satellite  oscillators  exceed  that  threshold. 

A  final  remark  briefly  challenges  a  noise  control  estimate.  The  estimate  is 
guided  by  all  three  analyses  (LIA),  (EA)  and  (SEA).  The  challenge  is  an  attempt  to 
dispel  an  assumption  that  often  leads  to  overestimation  of  the  benefits  that  may  be 
accrued  were  the  noise  control  guide  implemented.  Focusing  on  Eqs.  (1),  (2),  (9), 

(25)  and  (28)  one  may  define  the  ratio  So(y) »  of  ^^*0  kinetic  energy  stored  in 

the  coupled  to  the  kinetic  energy  Kk  (y)  stored  in  the  uncoupled  master  oscillator, 
in  the  form 

^oiy)  =  [e„kWE°okW]  =  [l  C(y)p  ] 

eo(.y)  =  P°(y){iM)inf(.y)l  •P‘’W  =  [n,(y)/n°M],  (39) 

where  /’‘^(jfis  the  ratio  of  the  external  input  powers  into  the  coupled  and  the 
uncoupled  master  oscillator.  The  loss  factor  rjf  (y)  is  stated  in  Eq.  (29);  namely 


rif(y)  =  [7oo(t)+7o5(t)];  Jios(y)  =  v^iy) 


(40) 
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The  blessings  of  the  noise  control  are  assessed  in  terms  of  the  smallness  of  as 

compared  with  unity;  a  beneficial  noise  control  is  one  for  which  s^iy)  «  1.  Often 

one  tends,  in  this  assessment,  to  argue  that  the  external  input  power  ratio  P^(y)  is 
largely  equal  to  unity  [7].  Using  Eq.  (40)  and  this  argument,  Eq.  (39)  reduces  to 

£o(y)  =>^o(y)  =  [Voo(y)/uHy)]  =  Vooiy) {^M+^ts(y)y^ ^  (4ia) 

Employing  the  first  order  approximation  of  ri%{y)  in  Eqs.  (40)  and  (41a),  one 
further  obtains 

£oiy)  =>  (t)  =  Vo  {Vo  +  ;  Piiy)  =  i-  (4ib) 

From  Eq.  (41)  it  becomes  clear  that  if  T]^{y)md  {rj^)  are  much  less  than  T]%iy) 
and  rfjiy),  respectively,  the  noise  control  that  can  be  projected  is  impressive.  For 

example,  as  Figs.  3-6  and  8-11  testify,  if  77o  =  (10“^)  [rjo^iy)  =  2(l+>’^)“’  •  (10“^).] 
then  for  strong  and  even  moderate  couplings,  the  projected  noise  control  is  two  to 
three  orders  of  magnitude.  In  passing,  except  for  a  minor  erosion  consideration, 
these  noise  control  benefits  are  independent  of  the  modal  overlap  parameter  (b), 
notwithstanding  that  the  undulations  at  low  values  of  (Z));  b  «\,  need  to  be 
recognized  and  handled  with  care.  An  obvious  question  then  arises:  Are  these 

projections  of  noise  control  real?  Is  the  ratio  P^{y)  or  its  first  order  approximation 

p/(y) ,  equal  to  unity  as  assumed?  After  all,  the  external  input  power  injected  into  a 
dynamic  system  is  critically  dependent  on  the  loss  factor  that  the  system  presents  to 
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the  external  drive.  As  Eq.  (29)  shows,  the  coupling  of  a  set  of  satellite  oscillators  to 
the  master  oscillator  enhances  the  loss  factor  that  the  master  oscillator  presents  to  the 

external  drive.  It  is,  thus,  expected  that  the  external  input  power  ratio  P'^iy)  harbors 
levels  that  exceed  unity,  thereby,  mitigating  the  noise  control  achievement  predicted 
in  Eq.  (41).  This  kind  of  mitigation  conforms  with  Le  Chatelier's  Principle  [18.]  If 
the  mass  of  the  satellite  oscillators  remain  within  reasonable  limits,  the  rosy  noise 
control  projections  are  mitigated  in  the  manner 

P°(y)  =  [I+HJW];  £^(y)  =  [1  +  H'(y)]«°(>')=  [/7co(>’)/'7jW]; 

=  (42a) 

and  the  equivalent  first  order  approximation 

^/(t)  =  [1  +  S7(t)];  <?o(T)  =  [l  +  S/(y)]£-/(T); 

=  (42b) 

where  use  is  made  of  Eqs.  (6),  (7),  (22),  (30)  and  (37)-(39)  [7,  17].  It  is  recalled  that 

El(y)  is  the  ratio  of  the  energy  stored  in  the  satellite  oscillator  and  in  the  couplings 
to  the  energy  stored  in  the  master  oscillator.  Design  conditions  may  be  found  to 
render  this  ratio  large  compared  with  unity  [7].  Except  for  minor  erosions  and 
undulations  that  may  need  to  be  recognized  and  handled  for  low  values  of  (b).  Figs. 
12-15  exemplify  that  for  strong  and  even  moderate  coupling  strengths  one  may  find 

levels  of  E^iy)  and  Hf(>’)  that  exceed  unity  by  an  order  or  two  of  magnitude  [7]. 
Such  large  levels  may  mollify  much  of  what  Eq.  (41)  promises.  Indeed,  Eqs.  (42a) 
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and  (42b)  suggest  that  if  rj^iy) ,  which  is  the  loss  factor  assigned  to  a  typical  satellite 
oscillator,  is  less  than  the  loss  factor  of  the  master  oscillator  in  isolation,  even 

a  noise  control  reversal  may  result  were  the  coupling  instituted.  The  unusual,  but 
possible,  conditions  for  such  a  reversal  are 

£oiy)  =  {vAy)lri"{y)}  >  h  Kiy)  »  l;  =  rj^(y)El{y)  »  (43a) 

£i(y)=^  s5(y)»i;  n^(y)^j(y)»^o,(y),  (43b) 

respectively  [7,  17].  [cf  Eqs.  (30)  and  (31).]  Of  course,  such  noise  control  reversal 
may,  when  it  occurs,  be  highly  surprising  to  those  who  factually  accept  the  validity 
of  Eq.  (41) . 
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Fig.  1.  Aset 


of  satellite  oscillators  coupled  to  a  master  oscillator. 


Fig.  2.  Loss  factors  j]{7)  assigned  to  the  satellite  oscillators  as  a  function  of  the 
normalized  index  (F)  and  with  the  mass  coupling  parameter  (m^)  equal  to  zero, 
a.  R  =  27. 


b.  R  =  7. 


Fig.  3.  Induced  loss  factor  T]l{y),  as  a  function  of  (y),  with  stiffiiess  control 
couplings.  [{MJN)  =  0.1  ;  modal  overlap  parameter  (b)  and  first  order  approx¬ 
imation  (FOA).] 

a.  (R  =  27)  Sprung-masses:  =1.0[a=0.0.],  g=mc=0. 

[Very  strong  coupling.] 


b.  (R  =  7)  Sprung-masses:  ac=1.0[a=0.0.],  g-  m^-0. 
[Very  strong  coupling.] 


Fig.  3.  Induced  loss  factor  rj^iy),  as  a  function  of  (y),  with  stif&iess  control 
couplings.  [{MJN)  =  0.1  ;  modal  overlap  parameter  (b)  and  first  order  approx¬ 
imation  (FOA).] 

c.  (R  =  27)  Satellite  oscillators;  =0.75  [a =0.25.],  g  = 

[Strong  coupling.] 

d.  (R  =  27)  Satellite  oscillators:  ac  =  0.15[Q:  =  0.85.],  g=  m^^O. 

[Moderate  coupling.] 
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Fig.  4.  Induced  loss  factor  rj^iy),  as  a  function  of  (y),  with  gyroscopic  control 
couplings.  [R  =  27  and  {Ms! M)  =  0.1 ;  modal  overlap  parameter  (b)  and  first  order 
approximation  (FOA).] 

a.  =  0  [<2=1.0.],  g  =  0.75 .  [Strong  coupling.] 

b.  ac  =  rhf;  =0  [<2=1.0. ],  g  =  0.15 .  [Moderate  coupling.] 


Fig.  5.  Induced  loss  factor  as  a  function  of  (y),  with  mass  control  couplings. 

[R  =  27  and  (M^/M)=0.1;  modal  overlap  parameter  (b)  and  first  order 
approximation  (FOA).] 

a..  ac  =  g  =0  [(2  =  1.5.],  fh^  =  0.15.  [Moderate  coupling.] 


h.  a^=  g  =0  [(2  =  1.03.],  -  0.03.  [Weak  coupling.] 
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a. 


Fig.  6.  Induced  loss  factor  ril{y),  as  a  function  of  (y),  with  mix  control  coupling 
forms.  [R  =  27  and  /M)  =  0.1 ;  modal  overlap  parameter  (b)  and  first  order 
approximations  (FOA).] 

a.  =  0.75  [« =1 .0.],  g  =  0.  [Very  strong  coupling.] 

b. a^=  0.53[«=.47.],  g  =  0.54,  m^=0.  [Strong coupling.] 


Fig.  6.  Induced  loss  factor  77j(y) ,  as  a  function  of  (y),  with  mix  control  coupling 
forms.  [R  =  27  and  {MJM)  =  0.1 ;  modal  overlap  parameter  (b)  and  first  order 
approximations  (FOA).] 

c.  =  0. 1  [a  =  .9.],  g  =  0. 1 1,  fhc  =  0.  [Moderate  coupling.] 

d. ac=  0.02[a=.98.],  g  =  0.22,  ^^=0.  [Weak coupling.] 


V\y) 


Fig.  7.  First  order  approximation  of  the  loss  factor  rf{y)  that  is  assigned  to 
individual  satellite  oscillators,  as  a  function  of  (y),  in  the  appropriate  frequency 
range.  [Modal  overlap  parameter  (b).]  [cf  Figs.  2a  and  2b.] 
a.  R  =  27. 


b.  R-7. 


(2.0)  (7;, 


Fig.  8.  Loss  factors  ri^{y)  (dashed)  and  r]l{y)  (solid),  as  functions  of  (y),  with 

stiffhess  control  couplings.  [(M^/M)  =  0.1  and  7o  =  (10“^);  modal  overlap  par¬ 
ameter  (b).]  [cf.  Fig.  3.] 

a.  [R  =  27]  Sprung-masses;  a^=1.0  [a  =  0.0.],  ^=mc=0.  [cf  Fig.  3a.]. 


b.  [R  =  7]  Sprung-masses;  a^=1.0  [a=0.0.],  g=mc=0.  [cfFig.  3b.]. 


d. 


Fig.  8.  Loss  factors  (dashed)  and  ril{y)  (solid),  as  functions  of  (y),  with 

stiffiiess  control  couplings.  [(M^ /M)  ==  0.1  and  7o=(10  );  modal  overlap  par¬ 
ameter  (b).]  [cf  Fig.  3.] 

c.  Coupling  as  specified  in  Fig.  3c.  [R  =  27.] 

d.  Coupling  as  specified  in  Fig.  3d.  [R  =  27.] 


Fig.  9.  Loss  factors  rj^iy)  (dashed)  and  rjHy)  (solid),  as  functions  of  (y)  with 

gyroscopic  control  couplings.  [R=275  (M^  / M)  =  1  and  77^  =  (10  ^) ;  modal  overlap 

parameter  (b).]  [cf  Fig.  4.]  (*) 

a.  Coupling  as  specified  in  Fig.  4a. 

b.  Coupling  as  specified  in  Fig.  4b. 


(*)  When  regions  of  curves  clearly  overlap,  the  color  of  the  one  with  the  higher 
modal  overlap  parameter  (b)  wins. 


Fig.  10.  Loss  factors  T]f{y)  (dashed)  and  (solid),  as  functions  of  (y),  with 

mass  control  couplings.  [R=27,  (M^ /M)  =  0.1  and  7^=  (10“^) ;  modal  overlap 


parameter  (b).]  [cf.  Fig.  5.]  (*) 

a.  Coupling  as  specified  in  Fig.  5a. 

b.  Coupling  as  specified  in  Fig.  5b 


Fig.  11.  Loss  factors  (y)  (dashed)  and  (solid),  as  functions  of  (y),  with 

mix  control  coupling  forms.  [R=27,  (M,/M)  =0.1  and  7^  =  (10“^) ;  modal  overlap 
parameter  (b).]  [cf  Fig.  6.]  (*) 

a.  Coupling  as  specified  in  Fig.  6a. 

b.  Coupling  as  specified  in  Fig.  6b. 
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Fig.  11.  Loss  factors  ?jfiy)  (dashed)  and  ri^y)  (solid),  as  functions  of  (y),  with 

mix  control  coupling  forms.  [R=27,  (M^  / M)  =0.1  and  7^  =  (10  ^) ;  modal  overlap 
parameter  (b).]  [cf  Fig.  6.]  (*) 

c.  Coupling  as  specified  in  Fig.  6c. 

d.  Coupling  as  specified  in  Fig.  6d. 


Fig.  12.  Modal  coupling  strength  a  function  of  (y),  for  stiffiiess  control 

couplings.  [{Ms/M)  =  0.\;  modal  overlap  parameter  (b)  and  first  order  approx¬ 
imation  (FOA).]  (*) 

a.  Coupling  and  number  of  satellite  oscillators  as  specified  in  Fig.  3a. 

b.  Coupling  and  number  of  satellite  oscillators  as  specified  in  Fig.  3b. 


Fig.  12.  Modal  coupling  strength  ^o(y) »  ^  function  of  (y),  for  stif&iess  control 

couplings.  [(M^/M)  =  0.1;  modal  overlap  parameter  (b)  and  first  order  approx¬ 
imation  (FO  A).]  (*) 

c.  Coupling  and  number  of  satellite  oscillators  as  specified  in  Fig.  3c. 

d.  Coupling  and  number  of  satellite  oscillators  as  specified  in  Fig.  3d. 


Fig.  13.  Modal  coupling  strength  Co(y) »  ^  function  of  (y),  for  gyroscopic  control 

couplings.  [R=27,  (M^/M)  =  0.1;  modal  overlap  parameter  (b)  and  first  order 
approximation  (FOA).]  [cf  Fig.  4.]  (*) 

a.  Coupling  as  specified  in  Fig.  4a. 

b.  Coupling  as  specified  in  Fig.  4b. 


a. 


b. 


Fig.  14.  Modal  coupling  strength  Cp(y)’  ^  function  of  (y),  for  mass  control 
couplings.  [R=27,  {M^l  M)  =  QA  \  modal  overlap  parameter  (b)  and  first  order 
approximation  (FOA).]  [cf  Fig.  5.]  (*) 

a.  Coupling  as  speeified  in  Fig.  5a. 

b.  Coupling  as  specified  in  Fig.  5b. 


Fig.  15.  Modal  coupling  strength  Co(y) ,  as  a  function  of  (y),  for  mix  control 
coupling  forms.  [R=27,  {MJM)  =  0.1 ;  modal  overlap  parameter  (b)  and  first  order 
approximation  (FOA).]  [cf  Fig.  6.]  (*) 

a.  Coupling  as  specified  in  Fig.  6a. 

b.  Coupling  as  specified  in  Fig.  6b. 


Fig.  15.  Modal  coupling  strength  as  a  function  of  (y),  for  mix  control 

coupling  forms.  [R=27,  (M^  /M)  =  0.1 ;  modal  overlap  parameter  (b)  and  first  order 
approximation  (FOA).]  [cf  Fig.  6.]  (*) 
c.  Coupling  as  specified  in  Fig.  6c. 


d.  Coupling  as  specified  in  Fig.  6d. 
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